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We discuss inverse problem of detection turbulence magnetic field helical properties using radio survey 
observations statistics. In this paper, we present principal solution which connects magnetic helicity and 
correlation between Faraday rotation measure and polarization degree of radio synchrotron emission. The 
effect of depolarization plays the main role in this problem and allows to detect magnetic helicity for certain 
frequency range of observable radio emission. We show that the proposed method is mainly sensitive to a 
large-scale magnetic field component. 
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1. INTRODUCTION 
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, Magnetic fields exist not only in compact astrophys- 
'ical objects such as planets and stars but also are ob- 
served everywhere in the Universe, interstellar space and 
can be attributed to galaxies and galactic clusters [T]. 
Substantially the dynamo theory explains the nature 
and evolution of cosmic magnetic fields 2 . The dynam- 
ics specific property of magnetohydrodynamics (MHD) 
systems are turbulent motions of media. The mean field 
'theory ^ predicts magnetic field generation, as a result 
of Qf-efFect which can appear under the condition of he- 
lical turbulent fiows. As a result magnetic field becomes 
helical too [Ij and can be described by magnetic helicity 
H = A-B, where A - is the vector potential of magnetic 
'field B = curl A. 

Recently the special role of magnetic helicity in space 
magnetic fields evolution processes is noted [5]. Total 
magnetic helicity of a system is integral of motion and 
conserved in the nondissipative limit. The results of 
theoretical and numerical researches show that the mag- 
netic helicity can accumulate in the system and suppress 
the generation mechanisms [6 . This put into question a 
L 5^ .possibility of the turbulent dynamo and has demanded 
construction of the adequate model describing dynamics 
H. 

The model of dynamo in the galactic disk has been 
added by equations describing outflux magnetic helic- 
ity that has allowed to overcome catastrophic suppres- 
sion dynamo processes [7 . Thus, mechanisms of solar 
dynamo have been considered and it is shown that al- 
lowance for the helicity of the small-scale magnetic fields 
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is of crucial importance in limiting the energy of the gen- 
erated large-scale magnetic field 8:. Using the results 
of [7 , it is proved necessity of coronal ejections for the 
strong large-scale solar magnetic field generation [9] . 

The development of the models that describe evolu- 
tion of magnetic helicity requires understanding of non- 
linear processes in multi-scale systems as well as notions 
about magnetic energy and helicity spectral distribu- 
tions, non-uniformity and anisotropy properties of the 
spatial distributions. It is extremely important to have 
factual material confirming presence of helicity and its 
connections with other components of the media. 

The observations of helicity in the solar convective 
zone indicate the existence of connection between the 
intensity of current helicity and dynamo processes jIO| . 
Study of MHD turbulence in laboratory conditions is 
extremely difficult (for review, see [IIH . Single success- 
ful experiments for measuring the turbulent magnetic 
fields il2j considerably differ in values of the character- 
istic parameters, primarily the magnetic Reynolds num- 
ber. Analysis of astrophysical observations remains the 
most promising direction of research in this question. 

In current astrophysical researches there is no gen- 
eral approach to derivation helicity of interstellar mag- 
netic fields. It is discussed a possible way to detect 
magnetic helicity from cosmic microwave background 
(CMB) fiuctuations data 13 . In some cases [T3], infor- 
mation about helicity can be extracted from the proper- 
ties of cosmic rays if their source is known. The authors 
of these researches noted that their approaches requires 
presence of high-accuracy observation data, which we 
don't have at present, therefore the practical applica- 
tion of these approaches is limited. 

New generation of radio telescopes (SKA and LO- 
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FAR) offers great opportunities [TS] because new high 
accuracy and resolution data about the space mag- 
netism win be available in the near future. Magnetic 
fields of the interstellar medium are the most suitable 
object for derivation MHD turbulence properties. Due 
to relatively large scales, we can neglect the contribution 
of regular magnetic fields, stars and planets and assume 
that the continuous electrically conductive interstellar 
medium is in a state of turbulent motion, which is ex- 
cited by explosions of supernovae [16J. Another signifi- 
cant feature of the interstellar medium is in its "trans- 
parency" , i.e. depth distribution of radio sources allows 
to make an analysis in all three dimensions. 

The aim of the paper is to show the possibility of 
magnetic helicity detection in the ionized interstellar 
plasma by statistical analysis of radio polarized obser- 
vations. We consider the model distribution of the mag- 
netic field with given properties and determine relation 
of magnetic helicity and correlation coefficient between 
Faraday rotation measure and polarization degree of ra- 
dio emission. 



2. INTERSTELLAR MEDIUM MODEL 

The simulation domain is a cube of side L. Let the 
coordinates x, y describe the sky plane and the axis z 
corresponds to the line of sight. For generation of artifi- 
cial polarized radio data distributions it is necessary to 
define some distributions of the ISM components such 
as magnetic field B, densities of relativistic ric and free 
thermal electrons n^. . An indicator of the magnetic field 
in the interstellar medium is synchrotron emission re- 
sulting from relativistic electrons passing through the 
magnetic field. 

At the first step of work we calculate three- 
dimensional distribution of uniform and isotropic mag- 
netic field. The input parameters of the model are 
power low of spectrum for magnetic energy distribution 
a, the turbulent energy scale I and magnetic helicity 
value. Also, the resulting magnetic field is solenoidal, 
i.e. V ■ B = 0. We can satisfy this conditions conve- 
niently using the Fourier representation of the magnetic 
field B expressed via the vector potential A 



B(k) = ik X A(k), A(k) = 



Ik X c| 



.|kr/2-i, (1) 



where k is the wave- vector, c = a+i b is the random 
complex vector, whose distribution determines magnetic 
helicity value. If the random vectors a and b have uni- 
form distribution the mean value of helicity (H) will be 
close to zero. If we choose only those pairs of vectors 



which give the same sign of k ■ (a x b) then (H) will 
be positive or negative, respectively. The extreme value 
of magnetic helicity for given magnetic energy will be 
obtained with 

, k X a . . 

The turbulent cells number along line of sight is de- 
fined by A'^ = [i^olj where ko — l^^ is the length of 
the wave-vector till which magnetic energy is equal to 
zero. Starting from ko, the energy spectrum obeys the 
Kolmogorov law a = —5/3. 

The next step of solving the problem is calculation of 
artificial polarized radio emission maps and RM images. 
The total intensity of synchrotron emission is given by 



Iix,y) = I e{x,y,z)dz, 

10 



(3) 



where e{x, y, z) is the synchrotron emissivity. Using sim- 
plified representations about spectral distribution of ric 
we can consider that e ~ nc{Bl + By). The synchrotron 
emissivity initially has some degree of polarization 7 
and the polarization angle is defined by a perpendicular 
direction to B in the sky plane {x,y). The intrinsic po- 
larization angle at the point of emission (x, y, z) is given 
by 

^o{x, y, z) = arctan(Sy/_B^) -I- 7r/2. (4) 

When polarized radio emission propagates thought mag- 
netized plasma, the polarization plane rotated by Fara- 
day effect. Thus, the polarization angle at some point 
with 



tpix, y, z) = iI;q{x, y, z) + A^RM(a;, y, z). 



(5) 



where A is the wavelength of observed radio emission, 
and RM is Faraday rotation measure which determined 
by the integral with variable upper limit 



RM(a:, y, 0) = X / neBz{x,y, z')dz' . 
Jo 



(6) 



Note that Faraday rotation of the polarization plane de- 
pends on the magnetic field component along the line 
of sight, while polarized intensity and polarization an- 
gle are defined by the perpendicular component. The 
observed Stocks parameters Q and U can be used to de- 
termine the complex intensity of the polarized emission 
P = Q + iU which given as 

P{x,y)=j e{x,y,z)exp{2ip{x,y,z)}dz. (7) 
Jo 

Superposition of the electromagnetic waves with differ- 
ent polarization angles causes depolarization, thus that 
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Fig.l. Radio maps. Distributions are represented in the observation plane {x,y): (a) Faraday rotation measure distri- 
bution 'R}A{x,y,z = L), (b) Polarization degree p for A = 0.05m and (c) for A = 0.2m. Image resolution is 256x256 
pixels. Black colour corresponds to minimum values and white corresponds to maximum values 



the observed polarization degree p = \P\/ 1 varies from 
to 7. Depolarization may be caused not only by phys- 
ical reasons and also by the limited radio telescope res- 
olution. In this work instrumental effects are not con- 
sidered. The physical size of the simulation domain is 
L = 0.5 kpc (1 kiloparsec « 3-10^^ ) that corresponds to 
the half-thickness of the galactic disk. The dimensional 
constant m ^ K = 0.81 if A is measured in meters, 
z in parsecs, and in cm^"^. We accepted the mean 
value of the magnetic field B = \ fiG and the thermal 
electron density = Icm^^ as a typical values for the 
ISM. 

Figure 1 shows typical view of calculated Fara- 
day rotation measures and polarization degree distribu- 
tions for wavelengths of radio emission 0.05 m, 0.2 m. 
These distributions of the radio data contain informa- 
tion about all magnetic field components. The changes 
of distribution details p for long A are explained by Fara- 
day depolarization. It is possible to see formation of 
thin black structures (figure Ic) in domains correspond- 
ing to the maximal values of Faraday rotation measure, 
this structures are typical for real astrophysical observa- 
tions. The analysis of these structures called "canals", 
allows to identify some properties of interstellar turbu- 
lence [T7]. Noted canals appearance refiects the fact of 
reliability for the chosen ISM model. 

3. STATISTICAL ANALYSIS 

Figure 2 shows probability distribution functions of 
RM for three levels of magnetic helicity. These distri- 
butions have symmetric shapes which are sufficiently 
approximated by the normal law. The difference be- 
tween distributions is insignificant that cannot be used 
for diagnostics of helicity. 
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Fig. 2. Probability distribution function of Faraday ro- 
tation measure for different levels of magnetic helicity 

The situation changes substantially if we consider 
joint probabilities. The joint probability distribution 
density of RM and p depending on sign of magnetic he- 
licity is shown in Fig. 3. Magnetic helicity destroys the 
symmetry of the distribution function. For (H) > 
(Fig. 3a) low degree of polarization is most likely cor- 
responds to negative values of RM, and for {H) < 
(Fig. 3c) - positive values. The quantitative estimation 
of revealed statistical characteristics can be derived us- 
ing the correlation coefficient 

^ ^ {RMp)-{RM){p) 

v/((i?A/2)_(i?M)2)((p2)_(p)2)' 

where mean value is taken in the observation plane (x, y) 
for z ^ L. The probability distribution function of C 
is defined by repeated calculations of random magnetic 
field realizations with given level of magnetic helicity. 
For construction of the results that shown in Fig. 4 we 
used 300 realizations for each level of magnetic helicity, 
respectively. The ranges of C values for each magnetic 
helicity level practically doesn't intersect that allows to 
identify the sign of magnetic helicity. As mentioned pre- 
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Fig. 3. Joint probability distribution of pair RM and p for three levels of magnetic helicity : (a) positive, (b) zero-order, 
(c) negative 
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Fig.4. Probability distribution function of C with A — 
0.2 m for three levels of magnetic helicity 
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Fig.5. Mean of correlation C depending on wave length 
A for different levels of magnetic helicity 

viously, the main reason for connection of RM and p 
is probably the depolarization effect caused by Faraday 
rotation, that confirms the relationship of mean corre- 
lation C with the wavelength A (see Fig. 5). 

For wavelengths less than 6 cm we have low Faraday 
depolarization that C is almost equal to zero. The most 



successful wavelength for given parameters of the inter- 
stellar medium model is A « 15 cm, which produces 
the extremum of C. The connection between C and 
the number of turbulent cells along the line of sight N 
also was investigated. Our calculations show that the 
value of C suddenly decreases with increasing number 
of turbulent cells (see Fig. 6). 

According to the mathematical modeling results, the 
magnetic helicity initiates a correlation between polar- 
ization degree and Faraday rotation measure. Thus 
strong correlation is approximately equal to 0.4 and 
achieved for a definite wavelength. Optimum wave- 
length for the observations will depend on the param- 
eters of the interstellar medium: the domain size L, 
magnetic field i?, the thermal electron density n^. How- 
ever, it is found that for maximum effect, the interstellar 
medium should provide specific Faraday rotation, rota- 
tion of the polarization plane through angle of 27r. And 
then relation 

RM A^ w L B ne A^ w 27r. (9) 

The correlation decreases with increasing of the num- 
ber of turbulent cells along the line of sight. It means 
that the proposed method for diagnosis of the magnetic 
helicity is mainly sensitive to large-scale magnetic field 
component. The noise appears during radio polarized 
observations process, probably, also produces accuracy 
problems for helicity detection. This influence can be 
assessed within the proposed model, but it makes sense 
to do so, when we have real data with known signal-to- 
noise ratio and other features. 

This work was supported by the Russian Founda- 
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and #07-01-96007-Ural) and the International Science 
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Fig.6. Mean of C as a function of turbulent cells num- 
ber N along the linc-of-sight for A = 0.2 m 
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